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SUMMARY 


The results of an experimental parametric investigation of whirl flut- 
ter are presented for a model consisting of a windmilling propeller-rotor, or 
"proprotor, ” having blades with offset flapping hinges mounted on a rigid pylon 
with flexibility in pitch and yaw. The investigation was motivated by the need 
to establish a large data base from which to assess the predictability of whirl 
flutter for a proprotor since some question has been raised as to whether flut- 
ter in the forward whirl mode could be predicted with confidence. In order to 
provide the necessary data base, the parametric study included variations in 
the pylon pitch and yaw stiffnesses, flapping hinge offset, and blade kinematic 
pitch-flap (63) coupling over a large range of advance ratios. Fifty cases of 
forward whirl flutter and twenty-six cases of backward whirl flutter are docu- 
mented. The measured whirl flutter characteristics, which include flutter 
speed, flutter frequency, direction of pylon whirl, and pylon yaw-to-pitch 
amplitude ratio and phase angle, are shown to be in good to excellent agreement 
with predictions from two different linear stability analyses which employ a 
two-dimensional, quasi-steady aerodynamic theory which neglects the effects of 
the unsteady wake, pitching moment, noncirculatory lift, and profile drag. On 
the basis of these results, it appears that proprotor whirl flutter, both for- 
ward and backward, can be predicted with linear stability analyses by using 
simple, two-dimensional, quasi-steady aerodynamics for the blade loading. 


INTRODUCTION 

The phenomenon known as propeller whirl flutter came under Intense investi- 
gation in i960 as a result of the loss of two Electra aircraft. Whirl flutter 
involves a self-sustained or divergent precessional motion or "whirl” of the 
propeller about its unperturbed position and can occur in a flexibly mounted 
propeller /nacelle installation in which the nacelle degrees of freedom allow 
the propeller plane to process in response to the gyroscopic moments associated 
with the rotating propeller. When precession occurs, aerodynamic forces and 
moments which can be destabilizing on the whirl motion are generated because of 
cyclic angle-of-attack changes on the blade elements of the propeller. Under 
some conditions the precession is in the same direction as the propeller rota- 
tion (forward whirl), and under other conditions the precession is opposite 
(backward whirl). On conventional (nonhinged) tractor propellers, it has been 
found that the instability invariably occurs in the backward whirl mode. 

During the period in which intense interest was being focused on whirl 
flutter for conventional propeller-driven aircraft, several VTOL aircraft which 
employed either articulated rotors or propellers having blades hinged to permit 
flapping were in research flight test programs. These aircraft employed the 
rotors or propellers for vertical flight and tilted the rotors or propellers 
forward approximately 90° for forward propulsion in the manner of an airplane. 
Because the rotors or propellers used on these aircraft perform the dual func- 



tlon of a lifting rotor and a propulsive propeller, the devices may be broadly 
referred to as propeller-rotors or "proprotors" for short. The terminology 
proprotor will be used in this report. Because of the flapping freedom of the 
blades, the whirl modes are affected, and a more complex variety of whirl flut- 
ter is possible. The concerns surrounding the Electra investigations provided 
the motivation to study the manner in which the whirl flutter characteristics 
of proprotors might be altered by the flapping freedom of the blades. Some 
early studies which ensued are reported in references 1 to 5. 

In reference 1 , the equations of motion for the whirl flutter analysis of 
a proprotor having offset flapping hinges were derived and employed in a sta- 
bility analysis of a 0.30-meter-diameter model which was used in a companion 
experimental study. For the case in which the blade flapping freedom was 
locked out, backward whirl was obtained and was in agreement with theory. For 
the flapping case, forward whirl was observed but the predicted instability was 
in the backward mode. Only by modifying the theory to include arbitrary phase 
lags in the blade aerodynamics in order to approximate the unsteady effects of 
the wake could a forward mode of instability be predicted. Reference 2 is an 
analytical investigation for a specific VTOL aircraft which utilized blades with 
offset flapping hinges . Reference 3 presented some experimental results for a 
different 0.30-meter-diameter model which had provision for two hinge offsets. 
Backward whirl was observed when the flapping was locked out and for the larger 
of the two hinge offsets. Forward whirl was observed for the smaller hinge off- 
set . In references 4 and 5 the theory of reference 1 was employed in a whirl 
flutter analysis of the model described in reference 3. The theory correctly 
predicted all the cases of backward whirl flutter but was unable to predict the 
forward whirl which was observed. Even by modifying the aerodynamic theory to 
include arbitrary phase lags to approximate the effects of the unsteady wake, 
forward whirl could not be predicted. 

While these investigations were being conducted, whirl flutter in a back- 
ward whirl mode was encountered during a full-scale wind-tunnel test of the 
Bell XV-3 convertiplane . The XV-3 had two large diameter, two-bladed proprotors 
which were mounted on tiltable pylons attached to the tips of the wings. The 
XV-3 had a teetering hub. As the name implies, the hub allows flapping of the 
blades, but the two blades flap about the shaft as a unit, one up, the other 
down, in the manner of a see-saw. References 6 eind 7 summarize the results of 
some analytical and experimental efforts aimed at explaining and correcting the 
whirl instability encountered on the XV-3. In particular, reference 6 indicated 
the possibility of both forward whirl and backward whirl for the XV-3, depending 
on the values of the system parameters. The possibility of both forward and 
backward whirl flutter was also indicated in reference 8 which examined the 
influence of flapping restraint on stability of a proprotor having centrally 
hinged blades by using a linear stability analysis and employing quasi-steady 
aerodynamics for the blade loading. However, no experimental results were shown 
to substantiate these predictions. 

Some more recent work directed at examining various aspects of the dynamic 
behavior of proprotor systems Is described in references 9 to 17. However, 
either because of the broad nature of the study or a lack of substantiating 
data, none of these studies were able to allay the skepticism which had developed 
regarding the ability to both predict forward whirl flutter and correlate thoge 
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predictions with experiment . With a view toward obtaining the experimental 
data needed for a more realistic assessment of the predictability of proprotor 
whirl flutter, a joint NASA/Grumman investigation of whirl flutter was conducted 
in the Langley transonic dynamics tunnel by employing an unpowered, 1 . 52 -meter- 
diameter model of a three-blade proprotor with offset flapping hinges and kine- 
matic pitch- flap (63) coupling. The proprotor was mounted on a rigid pylon 
which was restrained in pitch and yaw by springs. To provide a large sind varied 
data base, a range of pylon pitch and yaw stiffnesses, hinge offsets, and kine- 
matic pitch-flap (63) coupling angles were investigated over a wide range of 
windmilling advance ratios. Fifty cases of forward whirl flutter and twenty-six 
cases of backward whirl flutter were clearly identified. Two linear stability 
analyses available at the time of the experimental investigation, both employing 
simple two-dimensional, quasi-steady aerodynamics for the blade loading, were 
used in a companion whirl flutter analysis. Some preliminary results of limited 
scope pertaining to both the experimental and analytical aspects of this inves- 
tigation were previously reported in references I8 and 19 . The purpose of the 
present report is twofold: first, to provide a complete documentation of the 

experimental whirl flutter results which were given only limited treatment in 
references 18 and 19; second, to more fully substantiate the validity of analy- 
ses which employ simple two-dimensional, quasi-steady aerodynamics for predict- 
ing proprotor whirl flutter. 


SYMBOLS 


General Symbols 


Physical quantities in this report are given in both the International 
System of Units (SI) and U.S. Customary Units. All measurements and calcula- 
tions were made in U.S. Customary Units. 


Ap 


'n 


c 

D 

e 


pylon yaw-to-pitch amplitude ratio in whirl flutter mode 

n-1 

aerodynamic integrals, I D dri (n = 1, 2, 3 , 5 ) 

Ju 

Ul 

section-lift-curve slope per radian 

.T12 

aerodynamic integrals, I dri (n = 1, 2, 3) 

•'ui 

blade chord, m (ft) 
diameter of proprotor, m (ft) 

offset of flapping hinge from shaft center line, m (ft) 


o»^Y o measured uncoupled natural frequencies of pylon in pitch and yaw 

’ with proprotor nonrotating and blades locked to forward flap- 

ping stops, Hz 
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p] identity matrix 

Vif* 

Jp flutter advance ratio, _£ 

nD 

m blade mass per unit length, kg/m (slug/ft) 

N number of blades 

n proprotor rotational speed, Hz 

R radius of blade (measured from center line of rotation to 

blade tip), m (ft) 

r local blade radius (measured from center line of rotation), m (ft) 

t time , sec 

V free-stream velocity, m/sec (ft/sec) 

Vf free-stream velocity at flutter, m/sec (ft/sec) 

* -\I(SFW 

63 blade kinematic pitch-flap coupling angle, deg 

Cp»?y viscous damping of pylon relative to critical damping in pitch and 
yaw, respectively 

viscous damping relative to critical damping of blade flapping motion 
n nondimensional radial coordinate, r/R 

Tii,ri2 integration limits for blade aerodynamic integrals 

p air density, kg/m 3 (slug/ft^) 

4)F pylon yaw-to-pitch phase angle in whirl flutter mode, deg 

proprotor rotational speed, rad/sec 
OJF frequency of pylon whirl at flutter, rad/sec 

wp»0)y pylon natural frequencies in pitch and yaw, fp q and fy q, corrected 
to reflect zero coning, rad/sec ’ ’ 

(’) derivative with respect to time, d/dt 

1^ ]] square matrix 

4 


«R ) 



{ } 

column matrix 

[ ]-■' 

matrix inverse 

+ 

denotes forward whirl 

, 

denotes backward whirl 


Symbols for Appendix A 

a-ijb-i flapping degrees of freedom of proprotor disc in longitudinal and 
lateral directions^ rad (see fig, A1) 


C(j)y>C^2 

h-j ,h2 


hi ,h2 


Ib 


pylon viscous damping coefficients in pitch and yaw, N-m-sec/rad 
( Ib-f t-sec/rad ) 

distance between proprotor hub and pylon pitch and yaw axes, respec- 
tively, m (ft) 

distance between pylon center of mass and pylon pitch and yaw axes, 
respectively, m (ft) 


flapping inertia of blade. 



r^m dr, kg-m^ (slug/ft^) 


Ip ^yy, Ip ^22 moments of inertia of pylon in pitch and yaw about its cen- 

’ ’ ter of mass, kg-m^ (slug-ft^) 

Ip flapping inertia of proprotor, 5 Ip, kg-m^ (slug-ft^) 


Kh 


hub flapping spring. 



N-m/rad (Ib-ft/rad) 


K(|)y,K(|,^ 

kg 

Mp,My 


Mr 


Sb 


pylon stiffnesses in pitch and yaw, N-m/rad (Ib-ft/rad) 
blade flapping spring, N-m/rad (Ib-ft/rad) 

mass of pylon effective in pitch and yaw, respectively, kg (slugs) 



R 

m dr, kg (slugs) 


rm dr, kg-m (slug-ft) 


Sr 


= NSp, kg-m (slug-ft) 



s 

3o 

Y 


complex eigenvalue 

angle of built-in coning (precone), rad 


Lock number, 


pacR^ 


■•B 


= 1 kg-m^/sec^ (slug-ft^/sec^) 


X inflow ratio, V/QR 

<t»y><|>z pylon degrees of freedom in pitch and yaw (see fig. A1), rad 


Symbols for Appendix B 

a-| ,a 2 nondimensional distance between proprotor hub and pylon pitch and 
yaw axes, respectively, in proprotor radii 

Ag = A5 - 2eA^ + e^A^ 

Cq.,Cqp pylon viscous damping coefficients in pitch and yaw, N-m-sec/rad 
(Ib-ft-sec/rad) 


H inflow ratio, V/JJR 


ip»iy 


mass moments of inertia of pylon about pylon pitch and yaw axes, 
kg-m2 (slug-ft2) 


I1 


I 2 


I 3 



mr^ dr, kg-m^ (slug-ft^) 
mr(r - e) dr, kg-m^ (slug-ft^) 
m(r - e)2 dr, kg-m^ (slug-ft^) 


K 


M 


*^2 ’^3 ’^4 


kg-m^ (slug-ft^) 

N I m dr + Mjjub where is mass of hub, kg (slugs) 

•'0 

degrees of freedoBF (see fig. B1) 
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m(r - e) dr, kg-m (slug-ft) 



e = e/R 

X complex eigenvalue 

N>i , ^2 uncoupled natural frequencies of pylon/rotor combination in pitch and 
yaw, respectively, nondimensionalized by rotor rotational speed 

\>3 nonrotating, uncoupled flapping natural frequency of proprotor disc 

nondimensionalized by rotor rotational speed 

T nondimensional time, J2t 

Primes denote derivatives with respect to T. 


APPARATUS 


Wind Tunnel 


The experimental investigations were conducted in the Langley transonic 
dynamics tunnel which is a continuous flow, single return, variable pressure, 
slotted-throat tunnel having a test section 4.8? meters (16 feet) square with 
cropped corners. The control room and test section walls are provided with 
large windows for close viewing of the model. Although the runs were conducted 
in air under near atmospheric conditions at free-stream Mach numbers less than 
0.30, the tunnel is capable of operation at stagnation pressures from near 
vacuum to slightly above atmospheric and Mach numbers up to 1.2. Either air 
or Freon- 12 can be used as a test medium. 


Model 

The model (fig. 1) employed in the whirl flutter investigation was an 
adaptation of a 1/4. 5-scale industrial model of a specific tilt-proprotor air- 
craft design (ref. 17) which incorporated blades with offset flapping hinges 
and kinematic pitch-flap (63) coupling to reduce blade flapping. To obtain 
flutter at low tunnel speeds, the whirl flutter model employed a special-purpose 
reduced-stiffness pylon- to-wing-tip restraint mechanism. The resulting pylon- 
to-wing flexibility was such that the proprotor/pylon combination was effec- 
tively isolated from the wing and the model support structure so that the pre- 
dominant freedoms were blade flapping and pylon pitching and yawing. The whirl 
flutter model was unpowered (that is, windmilling) and incorporated features to 
permit variations of selected parameters to provide a wide range of whirl flut- 
ter configurations. Parameters which were adjustable include the pylon pitch 
and yaw stiffnesses, hinge offset, and blade kinematic pitch-flap (63) coupling. 

The range of the adjustable parameters tested represents a broad spectrum 
of practical values. Pylon pitch and yaw frequencies, nondimensionalized by 
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the proprotor rotational speed, varied from 0.23 to 2.04. The ratio of pylon 
pitch-to-yaw stiffness varied from 0.25 to 1.85. Hinge offsets of 0.05R and 
0.13R and kinematic pitch-flap (63) coupling angles of +6.75°, +10.5°, +20°, 
and +30° were tested, although not in all possible combinations. Emphasis was 
placed on the configuration having a hinge offset of 0.05R and a 63 of +20° 
because these values were nearest those which might be employed in practice. 
Values of 63 less than +6.75° could not be obtained because of mechanical 
interference problems and model design limitations. Structural damping of the 
pylon, which was not adjustable, varied from 2^p = 0.005 to 2?p = 0.032 in 
pitch and from 2^y = 0.021 to 2^y = 0.065 in yaw. The inflow ratio at flut- 
ter Vp/flR varied from 0.34 to 1.36. 

Design and construction .- The wing spar (fig. 2) consisted of a hollow 
aluminum beam with spanwise flanges which was attached to a steel plate at its 
inboard (root) end. Two wing struts extending from the tip of the wing spar to 
the lower edge of the steel plate provided a high level of wing vertical bend- 
ing stiffness. The steel plate was attached to a support structure consisting 
of a tripod arrangement of tubular steel beams which was bolted to the tunnel 
wall and floor. A plywood panel (fig. 1) passing through the vertical plane of 
symmetry of the model and attached to the tripod support structure was employed 
both to provide a seal for the open backside of the semi-fuselage and to serve 
as a reflection plane. For the whirl flutter investigation, the nonstructural, 
segmented aerodynamic panels which provided the spanwise and chordwise distribu- 
tion of wing airfoil contour were removed (fig. 2). This arrangement destroyed 
the wing aerodynamic lift and further decoupled the pylon motions from the wing. 
For convenience, however, the flaps (which were not germane to the whirl flutter 
investigation) were left on during testing of the configurations having a 
5-percent hinge offset. 

To obtain flutter at low tunnel speeds, the pylon was soft-mounted to the 
wing spar by means of a reduced-stiffness pylon-to-wing-tip restraint mechanism 
(fig. 3) which permitted independent variations in the pylon pitch and yaw stiff- 
nesses. The natural frequencies of the wing/pylon combination (with the pylon 
locked to the wing tip) in the fundamental wing beamwise, chordwise, and tor- 
sional modes were 29.6 Hz, 19.7 Hz, and 39.5 Hz, respectively. With the pylon 
unlocked and soft-mounted to the wing, most of the combinations of pylon pitch 
and yaw stiffnesses resulted in frequencies of the proprotor/pylon combina- 
tion which were well below 8 Hz and thus well removed from the lowest wing fre- 
quencies. It was verified, both experimentally and analytically, that even for 
the maximum pylon support stiffness, coupling between the proprotor/pylon system 
and the wing was negligible. Thus, the wing was effectively a rigid support 
structure for the proprotor/pylon system and the predominant freedoms were pitch- 
ing and yawing of the pylon and flapping of the blades. 

The pylon support mechanism (fig. 3) consisted of a boxlike housing which 
contained bearings allowing for freedom in pitch and yaw and two pairs of steel 
or aluminiun bars oriented to provide restraint in pitch and yaw. The arrange- 
ment of the bearings was such that the pitch and yaw axes were noncoincident. 

(The yaw axis was 0.043 m (0.1 4 ft) ahead of the pitch axis.) All combinations 
of the steel and aluminum bars used to provide the pylon pitching and yawing 
restraints were shown to have linear load-deflection (and hence stiffness) char- 
acteristics by loading and unloading them in a bench test. The hysteresis (and 
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hence damping) observed during the loading and unloading was always greater for 
the yaw spring than for the pitch spring. The pitch and yaw stiffnesses were 
independently adjustable by means of a slider which could be moved along the 
length of either the pitch or yaw spring and locked in position; thereby the 

effective length (and thus stiffness) of the spring was varied.- Damping in the 

pitch and yaw directions could not be set or varied independently but generally 
varied with changes in stiffness . 

The blades consisted of a negatively twisted (-23°) aluminum spar having 
spanwise distributed beamwise, chordwise, and torsional stiffnesses and covered 
with a continuous (nonsegmented) balsa wood skin bonded to ribs attached to the 
spar of the blade. Lead weights were bonded to the spar to obtain the proper 
weight and balance distribution. The desired level of blade Inplane stiffness 
was achieved by means of aluminum spar caps bonded to the upper and lower sides 
of the blade over the Inboard third of each blade and by several layers of boron 
filament tape bonded to the upper and lower surfaces of the blade from the root 
to about the 80-percent radius. This modification resulted in a blade first 
elastic nonrotating inplane frequency of about 28 Hz. Since the maximum rotor 

rotational speed in the whirl flutter tests was 17 Hz, the minimum blade inplane 

frequency on a per revolution basis was 1.65. The first beamwise cantilevered 
elastic mode frequency of the blade was 18 Hz, and the first cantilevered tor- 
sional elastic mode frequency was 66 Hz. 

A schematic illustration of the model hub and control system geometry 
is shown in figure 4. Each blade is attached to the hub assembly by means of 
a steel pin which forms the hinge about which the blade flaps. No flapping 
restraint is employed. The offset of the flapping hinge is varied by inserting 
spacers of different lengths between the inner part of the hub which is attached 
to the rotor shaft and the outer part which retains the flapping hinge. The 
blade pitch change axis is coincident with the centroidal axis of the spar which 
lies along the blade quarter-chord line. The pitch of the blades is varied col- 
lectively by means of a swashplate ahead of the hub which is linked to the pitch 
horns by pitch links. 

Blade kinematic pitch-flap coupling is provided by having the pitch link 
attachment point on the pitch horn offset radially from the flapping hinge (the 
distance d-| in fig. 4). For a leading-edge pitch horn arrangement, this geom- 
etry establishes a virtual flapping hinge so that when the blade flaps forward, 
the blade pitch decreases. The acute angle formed by the line defining the 
virtual flapping axis with the flapping hinge is the pitch-flap coupling angle 
63 . This angle is defined to be positive if the blade pitch decreases when the 
blade flaps forward, as in the present situation. The 63 angle is varied man- 
ually by moving the end of the pitch horn inboard (to reduce 63 ) or outboard 
(to increase 63 ). 

The model for the 13-percent hinge-offset configuration (figs. 5 and 6 ) 
differed somewhat from that of the 5 -percent hinge-offset configuration, (figs. 1 
and 3). These differences were due to the fact that the configuration having 
the 13 -percent hinge offset could be obtained only by increasing the radius of 
the proprotor. This requirement necessitated the use of a new rotating swash- 
plate and pitch link arrangement, a larger spacer, and several other minor 
changes in the control mechanism. To accommodate the increased blade radius, 
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it was necessary to remove the semi-fuselage. The wing flaps, which were left 
on during the testing of the 5-percent hinge-offset configuration, were also 
removed . 

Physical properties .- The physical properties of the model proprotor and 
pylon for both hinge offsets are summarized in table I. The pylon inertial 
properties in pitch and yaw for each of the hinge offsets are different because 
of the noncoincidence of the pylon pitch and yaw axes. The rotating swashplate 
and pitch links move back and forth on the pylon as collective pitch is changed 
but affect the pylon inertias about the pitch and yaw axes by less than 1 per- 
cent. This variation was neglected in the analyses and is not reflected in val- 
ues given in table I. 

The pylon support stiffnesses in pitch and yaw, which were the parameters 
most varied in the test program, were determined by "plucking" the model in 
pitch and yaw. Prior to each run, the blades were positioned against, and 
secured to, the forward flapping stops at a 13.5° coning angle and the model 
was plucked first in pitch and then in yaw and the subsequent decaying oscilla- 
tion in each degree of freedom recorded on a direct-write oscillograph. The 
frequencies of the resultant transient time histories were then used in combi- 
nation with the known inertial properties of the model to get the pitch and 
yaw spring rates. The damping, as a fraction of critical, was also obtained 
from the traces of the decaying oscillations by using the familiar logarithmic 
decrement method. Typical traces of decaying oscillations obtained by plucking 
are shown in figure 7. The plucks were repeated after each run in order to 
identify any large change or structural deterioration that may have occurred 
during the run. In general, only negligible changes in both frequency and damp- 
ing were identified by these comparisons. The measured values of frequency and 
damping shown in table II represent an average of the pre- and post-run values. 

Instrumentation . - Blade flapping relative to the hub plane was measured on 
one blade by means of a soft flexure instrumented with strain gages which was 
clamped to the root of the blade and reached across the flapping hinge. Blade 
beamwise, chordwise, and torsional moments were measured on one blade by means 
of strain gages which were mounted on the spar of the blade at the root and at 
the midspan position. Data measured in the rotating system were brought out 
through the transmission to the fixed system by means of a slip-ring assembly. 
Pylon motions in pitch and yaw were measured by two accelerometers located on 
the forward part of the pylon and well removed from the axes of pitch and yaw. 
One accelerometer was oriented to sense pitch, and the other accelerometer was 
oriented to sense yaw. Rotor rotational speed was indicated directly on a 
tachometer driven by the signal from a magnetic pickup on the model. The azi- 
muthal position of an arbitrarily selected reference blade was given by the 
same pickup which pulsed every 90° of rotor revolution. The pulse correspond- 
ing to a specific azimuth position of the reference blade was larger than the 
other three. This feature was used mainly for dynamically balancing the rotor. 
Electrical controls were used to vary blade collective pitch (and hence rotor 
rotational speed) remotely from the tunnel control room. Data were recorded 
on an oscillograph and on magnetic tape. A short motion picture was taken of 
each flutter point. 
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TEST PROCEDURE 


For a given hinge offset and pitch-flap coupling angle, the pylon support 
stiffnesses and damping in pitch and yaw were first determined by plucking the 
pylon as described above. Whirl flutter points were then experimentally estab- 
lished in the following manner: A preselected windmilling rotor rotational 

speed was set at some low tunnel speed by remotely adjusting the collective 
pitch of the blades. The model was then transiently excited by plucking a light- 
weight cable which was attached to the pylon (fig. 3 ) and routed to the control 
room. The stability of the response was observed both visually and on an oscil- 
lograph. If the model was stable, the tunnel speed was increased by a small 
amount, the collective pitch retrimmed to produce the same rotor speed, and the 
model again excited. This procedure was repeated until a sustained, approxi- 
mately constant amplitude whirl oscillation (indicating neutral stability) was 
observed. The tunnel conditions were then recorded as the conditions for flut- 
ter and a short motion picture of the whirl instability taken. Tunnel speed 
was then reduced, the rotor speed was set to a new value, and the procedure 
was repeated. By proceeding in this way, several whirl flutter points could 
be obtained for a single setting of the pylon support stiffnesses but at vari- 
ous values of rotational speed. 

The exact point of neutral stability was established in most cases. In 
those cases where neutral stability was not achieved, the flutter point was 
established from the last few decaying (or diverging) traces by extrapolation. 

The cable used to excite the model was also used to restrain the model in those 
instances in which the pylon motions were divergent. 


DESCRIPTION OF ANALYSES 

Two linear stability analyses based on two different mathematical models 
(refs. 1 and 15) were employed in the whirl flutter analysis of the model con- 
figurations tested. Both analyses assume that the proprotor has three or more 
blades and is mounted on a rigid pylon which is restrained elastically in pitch 
and yaw by springs. The proprotor is assumed to be windmilling (unpowered) in 
an axial-flow condition and to consist of three or more rigid blades which are 
restrained elastically in flapping by springs. The blades are assumed to have 
constant chord and to incorporate kinematic pitch-flap ( 63 ) coupling. The dynam- 
ics of the present model are taken to be expressed in terms of only four of the 
eight degrees of freedom considered in references 1 and 15: pitch and yaw of 

the pylon and two cyclic flapping modes. To an observer in the nonrotating sys- 
tem, the two cyclic flapping modes are patterned so that the tip path plane 
appears to be pitching and yawing; that is, the rotor appears to be flapping 
longitudinally and laterally. For a proprotor having three or more blades and 
operating in axial flow, the use of these tip-path-plane coordinates as degrees 
of freedom in lieu of the individual blade flapping coordinates as degrees of 
freedom is an expedient which allows for the removal of the periodicity from 
the equations of motion. The introduction of these coordinates amounts to a 
coordinate transformation and is discussed in both references 1 and 15. 

A two-dimensional, quasi-steady aerodynamic theory which neglects the 
unsteady wake, pitching moment, noncirculatory lift, and profile drag is 
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employed to calculate the blade aerodynamic loading. The quasi-steady lift 
is taken to be completely determined by the resultant of the velocity components 
acting tangential and perpendicular to the blade chord at the quarter-chord of 
the section. The flow is assumed to be incompressible and the section-lift- 
curve slope to be constant over the length of the blade. 

The first analysis which was used to predict the whirl flutter behavior 
of the model is based on ttet developed in reference 15 for a gimbaled propro- 
tor. A gimbaled proprotor is characterized by blades which are rigidly attached 
to a hub assembly which is connected to the shaft by a gimbal (universal joint) 
housed within the hub. The blades of a gimbaled proprotor flap about the shaft 
by virtue of the gimbal mounting arrangement. Since the model tested had an 
offset flapping hinge, the inertial and aerodynamic restoring moment effects 
associated with the offset hinge were approximately accounted for in the gim- 
baled analysis by determining an equivalent hub spring which preserves the 
rotating flapping natural frequency of the hinged blade with aerodynamic forces 
included. The equations of motion used in the whirl flutter analysis are sum- 
marized in appendix A. The second analysis which was applied to the model is 
based on the derivation given in reference 1 which was developed explicitly for 
a proprotor having blades with offset flapping hinges. The equations of motion 
with four degrees of freedom appropriate to the present case are given in appen- 
dix 5 of reference 1 as a special case of the more general development. Several 
additions to those equations are made herein, however, in order to make them 
applicable to the present model. These additions include pylon inertia and 
damping about the pitch and yaw axes, flap-hinge damping, and a large pitch-flap 
coupling angle. The limits of integration for the blade aerodynamic integrals 
are also changed from those given in reference 1 . The resultant equations are 
summarized in appendix B. 

The steady-state coning angle of the blades during wind-on operation with 
the rotor in a windmilling condition was always negative (downwind) but never 
exceeded 1°. By use of the gimbaled proprotor analysis which accounts for blade 
built-in coning (precone), it was established analytically that these small con- 
ing angles have a negligible effect on the predicted whirl instability. Thus, 
for convenience in the adaptation of the theory of reference 1 (which did not 
include steady-state coning) to the present investigation, all the whirl flutter 
analyses were made by assuming the steady-state coning to be zero. The more gen- 
eral version of the gimbaled proprotor analysis which includes the wing degrees 
of freedom and is described in reference 15 was used to analytically verify that 
as a consequence of the soft pylon mounting arrangement in combination with a 
stiff wing support structure, the degrees of freedom of the wing had a negligi- 
ble effect on the predicted whirl flutter characteristics of the model. 

The structural damping about the blade flapping hinge was not measured and 
was taken to be zero in the analyses. This assumption was judged to be reason- 
able because of the low friction at the flapping hinges. Since the tunnel den- 
sity varied slightly from run to run, an average value of 1.23 kg/m3 (0.00238 
slug/ft3) was used in the analyses. The blade lift-curve slope was assumed to 
be constant and equal to 5.73 for the analytical studies. The blade aerodynamic 
integrals were evaluated by using a lower limit of integration which reflects 
the fact that because of root cutout, the lifting portion of the blade begins 
8.6 cm (3.^ in.) outboard of the hinge location, and by using an upper limit 
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which reflects the assumption that a tip length equal to one-half of the chord 
develops no lift. These considerations lead to limits of integration of 0.1 6 R 
and 0.9^R for the 5-percent hinge-offset configurations and limits of 0.24R and 
0.94R for the 13-percent hinge-offset configurations. 

In the computer implementation of the analyses, the measured natural fre- 
quencies of the model in pitch and yaw, fp q and fy o> ss part of 

the input data and corrected to reflect a zero coning condition. The resulting 
frequencies, in nondimensional form, are given by Up/fl and Wy/fl in table III. 
The blade kinematic pitch-flap coupling angle 63 varies slightly with changes 
in collective pitch by virtue of the varying geometry which defines this angle. 
These small changes are accounted for in the computer programs. The equivalent 
hub spring which is employed in the glmbaled proprotor analysis to account for 
the effects of the offset flapping hinge is a function of the blade physical 
properties, hinge offset, airspeed, and rotational speed. Since the flutter 
points identified with the model represent a large number of V,J2 combinations, 
the equivalent hub spring was calculated in the computer program. 


RESULTS AND DISCUSSION 
Comparison of Theory and Experiment 

A summary of all the model configurations which exhibited a whirl flutter 
instability along with the pertinent data defining each flutter condition is 
given in table II. A companion summary plot of all the experimental flutter 
points identified with the model is given in figure 8 where the data are plotted 
as a function of the pylon yaw and pylon pitch uncoupled natural frequencies 
corrected to zero coning and nondimensionalized by the rotor rotational speed. 
The data are seen to be grouped into three distinct regions emanating from the 
origin: a central region associated with symmetric and nearly symmetric pylon 

frequencies, and two side regions associated with nonsymmetric pylon frequencies 
which have a boundary with the center region at pylon pitch-to-yaw frequency 
ratios of about 0.75 and 1.5. For configurations having pylon frequency ratios 
falling in these nonsymmetric regions, flutter was always in the backward whirl 
mode. All the cases of forward whirl flutter fall near the line of equal pitch 
and yaw frequencies. Ten cases of backward whirl flutter also fall in the cen- 
tral region. Eight of these points represent all the whirl flutter points iden- 
tified with the configuration having a 13-percent hinge offset. The other two 
points are the two cases of backward whirl identified in run 48 for a 5 -percent 
hinge offset. 

A summary of the analytical results that correspond to the experimental 
results shown in table II is given in table III. The flutter phase angle (|)f 
which is obtained from the solution of the equations of appendix A is formally 
different from that obtained from a solution of the equations of appendix B. 

This difference is due to the difference in definition of the positive sense 
for the pylon pitch degree of freedom. (Compare figs. A1 and B1 . ) The phase 
angles from the two theories are related according to 

<l>Flref. 1 = <t>Flref. 15 + 1®°° 
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To provide for a direct comparison of these phase angles, the phase angles given 
in table III for the theory of reference 15 have been transformed to the nota- 
tion of reference 1 by using the preceding relation. 

A run-by-run comparison of both theories with experiment using the informa- 
tion provided in tables II and III is made in figure 9 for the flutter speeds 
and in figure 10 for the associated flutter frequencies. The flutter speeds 
shown are nondimensionalized by the tip speed of the rotor to give the flutter 
inflow ratio Vp/flR. This ratio is related to the flutter advance ratio Jp 
according to Jp = irVp/S2R. The flutter frequencies are nondimensionalized by 
the rotor rotational speed to give the frequencies on a per-revolution basis. 
Here, and in subsequent figures, the flutter results are given as a function of 
the pylon pitch uncoupled natural frequency corrected to zero coning and nondi- 
mensionalized by the rotor rotational speed. For the configurations of runs 48, 
50, and 65, the analyses indicated that there were two flutter boundaries which 
were close to the measured flutter boundary. Both of these analytical bounda- 
ries are shown in those cases. Both theories are in good to excellent agreement 
with experiment for both flutter speed and flutter frequency. The direction of 
pylon whirl at flutter was established both visually during the test and by 
inspection of the pylon pitch and yaw traces on the oscillograph records. The 
whirl direction predicted by theory was in every case in agreement with that 
observed experimentally. 

Flutter in the backward whirl mode occurred at low frequencies and in what 
might be termed a "rotor mode" since the flutter frequencies were near those of 
the low-frequency cyclic flapping mode with the pylon degrees of freedom locked 
out. Flutter in the forward whirl mode generally occurred at frequencies which 
were higher than those associated with flutter in the backward whirl mode and in 
what might be termed a "pylon mode" since the flutter frequencies were near the 
natural frequencies of the system in pitch and yaw with the flapping degree of 
freedom locked out. These aspects of the whirl dynamics of proprotor /pylon sys- 
tems are discussed in references 6, 7, and 15. The flapping in a pylon mode is 
generally large relative to the shaft (but small relative to space) whereas the 
flapping is a rotor mode is generally small relative to the shaft (but large 
relative to space). ^ Typical mode shapes which illustrate this flapping behav- 
ior for a proprotor/pylon system with four degrees of freedom and having hinged 
blades are shown in references 1 and 4. Both the backward whirl and forward 
whirl instabilities were accompeuiied by divergent motions in several instances 
in which the flutter speed was inadvertently exceeded. There were also several 
Instances in which both the backward and forward whirl instabilities were pre- 
cipitated merely by Increasing the tunnel speed, no external excitation via the 
lightweight plucking cable being required. 

A comparison of the measured pylon yaw-to-pitch amplitude ratios and phase 
angles at flutter with those predicted using the theory of reference 1 is given 
in figures 11 and 12. Reading accuracy of the phase angle from the oscillograph 
records was about ±10°. The phase angles in the whirl mode were not 90° (as 
they would be in a vacuum) but varied between about 60° and 120° because phase 
differences other than 90° are induced on the pylon whirl motion by the aerody- 
-** - . . - - 

^The XV-3 instability was a backward whirl associated with the low- 
frequency rotor flapping mode. 
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namic loading acting on the blades if the model pylon is not symmetric in the 
pitch and yaw directions. Similar comparisons, not shown, were obtained by 
using the theory of reference 15. Thus, in addition to correctly predicting 
the whirl directions, the theories of references 1 and 15 also predicted the 
amplitude ratios and phase angles reasonably well. 


Effect of Design Parameters 

Some selected composite results which show the effects of pitch-flap cou- 
pling, hinge offset, and asymmetry of the pylon support frequencies (stiffnesses) 
on whirl flutter stability are shown in figures 13 to 15. The effect of varia- 
tion of 63 on stability is given in figure 13 for a configuration having nomi- 
nally symmetric pylon pitch and yaw frequencies (within 5 percent of each other) 
and 5-percent hinge offset. The results show a strong decrease in flutter 
Inflow ratio (and hence flutter speed for a fixed rotor rotational speed) with 
increasing positive 63 . Although only positive values of 63 were investi- 
gated, reference 15 has shown that negative values of 63 are equally destabi- 
lizing on whirl flutter. Flutter is in the forward whirl mode except for the 
two points, denoted by solid symbols, which are in the backward whirl mode. The 
analytical results shown are based on the assumption of a symmetric pylon fre- 
quency configuration which is the mean of the pylon frequencies of the data 
points included in figure 13. Similarly, the damping used is the mean of the 
pitch and yaw values. 

The effect of hinge offset on stability is illustrated in figure 14 for 
configurations having nominally symmetric pylon pitch and yaw frequencies and 
63 = 20°. Increased hinge offset is seen to be strongly stabilizing. Theoreti- 
cal results are again based on the use of mean values for both the pylon support 
frequencies and dampings. 

The effect of asymmetry in the pylon support frequencies is indicated in 
figure 15 for the configuration having 5-percent hinge offset and 63 = 20 °. 
Again, the symmetric data represent nominally symmetric configurations in which 
the pylon frequencies are within 5 percent of each other whereas the nonsymme- 
tric data include configurations which satisfy the indicated frequency asymme- 
tries. The flutter inflow ratio is plotted against the nondimensionalized pylon 
pitch natural frequency except for the cases in which Wy < o)p, where the non- 
dimensionalized yaw frequency is shown for convenience in plotting. The results 
confirm earlier results (refs. 9 and 15) that for highly unsymmetric configura- 
tions in the pylon support frequencies (and hence stiffnesses) , whirl flutter 
occurs in the mode corresponding to the lower pylon frequency (stiffness) and 
increasing the higher frequency (stiffness) does not increase the flutter speed. 
Flutter is in the backward whirl mode for all the asymmetric conditions. 


Analytical Damping Studies 

It was stated earlier that the damping at the flapping hinge was not mea- 
sured and was taken to be zero in the analyses of the model configurations which 
exhibited whirl flutter. Reference 16 pointed out that internal damping, that 
is, damping which is exhibited in the rotating system, can be destabilizing on 
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proprotor/pylon whirl stability. It was also stated earlier that the pylon 
damping in pitch and yaw was not adjustable and could not be varied indepen- 
dently of the pitch and yaw stiffnesses. In order to analytically identify 
the effect of hinge damping and pylon damping on whirl flutter stability of 
the present model, some analytical studies in which these parameters were 
varied were made on three of the model configurations. The results of these 
studies are given in figures 16 to 21. 

The effects of hinge damping and pylon damping on whirl stability of a 
symmetric pylon frequency configuration in which instability occurred in the 
forward whirl mode are shown in figures 16 and 17, respectively. From fig- 
ure 16, hinge damping is seen to reduce the inflow ratio at which the forward 
whirl mode becomes unstable, the reduction becoming more significant with 
decreasing pylon support frequency. A reduction in pylon damping (fig. 17) 
has a moderate destabilizing effect on stability. 

The effects of hinge damping and pylon damping on whirl stability of an 
unsymmetric configuration in which instability occurred in a backward whirl 
mode are shown in figures 18 and 19, respectively. Neither hinge damping nor 
pylon damping is seen to have any significant effect on the backward mode of 
instability associated with an uns3rmmetric pylon frequency configuration. 

The effects of hinge damping and pylon damping on whirl stability of a 
symmetric frequency configuration in which instability occurred in a backward 
whirl mode are shown in figures 20 and 21. It is seen that, again, neither 
hinge damping nor pylon damping has any significant effect on . the backward mode 
of instability associated with a symmetric pylon configuration. 


CONCLUSIONS 

Results of an experimental parametric investigation of proprotor/pylon 
whirl flutter conducted in the Langley transonic dynamics tunnel are presented 
for a model consisting of a proprotor having blades with offset flapping hinges 
mounted on a rigid pylon with flexibility in pitch and yaw. The parametric 
study included variations in the pylon pitch and yaw stiffnesses, offset of 
the flapping hinge, and blade kinematic pitch-flap (62) coupling over a large 
range of windmilling advance ratios. Fifty cases of forward whirl flutter and 
twenty-six cases of backward whirl flutter are documented. Measured whirl flut- 
ter characteristics, which include flutter speed, flutter frequency, direction 
of pylon whirl, and pylon yaw-to-pitch amplitude ratio and phase angle, are 
shown to be in good agreement with predictions from two different linear stabil- 
ity analyses which employ two-dimensional, quasi-steady aerodynamics in which 
the effects of the unsteady wake, pitching moment, noncirculatory lift, and 
blade profile drag are neglected. On the basis of the results shown, the fol- 
lowing conclusions are indicated: 

(1) Proprotor whirl flutter, both forward and backward, can be predicted 
with linear stability analyses which employ simple, two-dimensional, quasi- 
steady aerodynamics for the blade loading. 
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(2) Stability analyses developed for a gimbaled proprotor can be applied 
to proprotors having blades with offset flapping hinges provided the blade 
inertial and aerodynamic terms associated with the offset of the hinge are 
accounted for through an equivalent hub spring which preserves the rotating 
flapping natural frequency of the hinged blade with aerodynamic forces included. 


Langley Research Center 

National Aeronautics and Space Administration 
Hampton, VA 23665 
October 12, 1977 
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APPENDIX A 


WHIRL FLUTTER EQUATIONS OF MOTION FOR A GIMBALED PROPROTOR 


The linear equations of motion for the whirl flutter analysis of a gimbaled 
proprotor are given in reference 15 in the general matrix form 

W{q} + [r]{q) + [c]{q} + [K]{q) = [AD]{q) + [AS]{q} (A1) 



ordering the four degrees of freedom shown in figure A1 
according to 


in the column matrix 




9y 




the square matrices appearing in equation (A1) are given b 


(A2) 


IR 


0 


[m] 


0 

^R + ^ SrPo^I 


IR 

0 


Ir + 1 Sp3ohi 0 

0 i ^R^oh2 

Mph-|^ + Ip + Mph-]^ 0 

+ 2Sp3ohl + Ip^yy 


Ir + I Sp3oh2 0 Myiip^ + Ip^2z + Ir 

+ Mph2^ + 2Sp3oh2 J 
(A3)2 


^The matrix expression given in equation (A3) is a special case of the 
more general eight-degree-of- freedom expression for Cm 3 given in reference 15. 
Six of the precone terms which appear in this expression are incorrect (low by 
a factor of 2). The result given by equation (A3) reflects this correction. 
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[ad] □ 


[r] 


APPENDIX A 


0 

2Ir 

0 

2Ir 

-2Ir 

0 

-2Ir 

0 

0 

2Ir 

0 

2Ir 

-2Ir 

0 

-2Ir 

0 


(A4) 


Ch 


[c] 




c<t> 


y 



(A5) 


Kh JiCn 

J2Ch Kh 

K<t>y 

K<t>z 


-A 5 /S 2 

AhX3q/^ 

-A^/^ - A2X2hi3o/J^R 

-A 3 Xh 2 /J^R + A4X3q/^ 

-A4X3q/S2 

-A^/J^ 

A 3 Xhi/fiR - AuX&q/SI 

-A^/Ji - A2X^h23o^^^ 

— A^/S 2 

Ai|X3o/i^ 

+ A 3 hiX/S 2 R 

-A5/n - A2X2h-|3o/^R 

— A-|X2hi2/J2R2 

— A3Xh2/J2R + Ai|X3o/^ 
+ A 3 hiX/f 2 R 

Ai^X^q/Q 

A‘^h2X/QY{ 

-A 5 /S 2 

A 3 Xhi/S 2 R - A2|X3o/S^ 
- A 3 h 2 X/QR 

-A^/n - A 2 X^h 26 o/SR 

- AiX 2 h 2 ^/S 2 R 2 


(A7) 
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-B3 tan ^3 - A4XBQ 

-A5 + XBqB 2 tan 63 

A2X3 Bo 

A3X2 


+A5 - XBqB 2 tan 63 

-B3 tan 63 - X3 qA4 

-A3X2 

A2X3 Bo 

[as] = 1C 

-B3 tan ^3 - XB0A4 
- A3hiX/R 

^®o®2 ^3 “ ^5 

+ (X/R)Bih-) tan 63 

A2^3$q 

+ AiX3hi/R 

A3X2 


A5 - XBqBp tan 63 
-(X/R)B-|h2 tan 63 

-B3 tan §3 - XB0A4 
- A3h2X/R 

-A3X2 

A2X3 $o 
+ AiX3h2/R 


(A8) 


To facilitate comparison of the present equations with those of reference 15, 
the notation of that reference is maintained as much as possible. 

The viscous damping coefficients Cfj, C,j) , and appearing in [cj and 

[k] are evaluated from the expressions ^ 


Cr = 2IrC3J^ \ h + 
C<t>y = 2CpUp,oIp,o 


(A9) 




where ^3 , Cp, and Cy are the damping ratios in flap, pitch, and yaw, respec- 
tively, Wp o = PTTfp^Q, ojy^o = PTTfy^o, and Ip^Q and ly^o are the mass 
moments of inertia of the model in pitch and yaw with the blades locked to the 
forward flapping stops., The equivalent hub spring KH,eq obtained from 
considerations given in appendix C. 

It should be noted that the hub damping coefficient Cr appears in the 
stiffness matrix [X] as well as in the damping matrix [cj . This is a conse- 
quence of the fact that the viscous structural damping associated with the flap- 
ping motion is regarded as being in the rotating system and the equations of 
motion are written with respect to a fixed or nonrotating system. 

Transferring the aerodynamic matrices in equation (A1) to the left-hand 
side and combining matrices yields the resulting equations of motion in the 
form : . 


[M]{q> . [C]{q> . [K]{q) = (o) 


(A10) 
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where 


[m] = [m] 

[c] = [r] - [ad] .[c] (AID 

[k] = [k] - [as] 


The matrix is symmetric and positive definite; [c] and [k] are nonsymmetric . 
By premultiplying equation (A10) through by 

[I]{q) * [M]-'[C](.^ . M-'[K](q> = (o) (S12) 

and introducing the four generalized velocities as auxiliary variables by means 
of the matrix identity 


{q> - {q> = {0> 

equations (A12) and (A13) can be combined to yield 



By defining 


0 

-M-'W 


■ [I] 

-W-'[c] 




equation (A14) can be written in the compact form 

{x>= CA]{x} 


(A13) 


(A14) 


(A15) 


(A16) 


where the definition of Ca] follows from equation (A14). Note that the four 
original second-order equations have been transformed to eight first-order equa- 
tions. By assuming a solution of the form 

{x} = {Xo}eSt (A17) 

equation (A16) reduces to 


[a](Xo]> = s(Xo} (A18) 

Equation (Al8) is in standard eigenvalue form and is amenable to solution by 
using standard eigenvalue techniques. 


The solution of equation (Al8) leads to eight complex eigenvalues Sp and 

eigenvectors ^Xq(p[}. Since [^A] is real, the eigenvalues and eigenvectors occur 
in four complex conjugate pairs. Only roots having positive imaginary parts 
need be considered since the roots having negative imaginary parts do not pro- 
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vide any new linearly independent solutions. The pth complex eigenvalue having 
a positive imaginary part has the general form 

Sp = cip + i3p ( A 1 9 ) 

and can be interpreted as 

Sp = -CpWp + - Cp^ (A20) 

where Cp is the damping as a fraction of critical damping and ojp is the 
undamped coupled frequency. For u)p ^ 0 if ?pfJp > 0> the motion is exponen- 
tially divergent. By plotting (?p,03p) as a function of velocity in the complex 
plane, a concise picture of the variation of system stability with airspeed can 
be established. 


The complex vector associated with the eigenvalue Sp can be written in 
the form 



(A21) 


The upper four elements in each modal vector of eight elements define 

the mode shape. Since these elements are complex, phase differences exist 
between the component harmonic motions in a given mode. The relative amplitude 
and phasing between the pylon pitching and yawing motions in a given mode can 
be ascertained by converting the complex elements of the mode to polar form and 
then normalizing on one of them. 
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Degrees of freedom 

a,^ Longitudinal flapping of proprotor disc 
Lateral flapping of proprotor disc 
0y Pylon pitch rotation 
(j>^ Pylon yaw rotation 



Figure A1.- Four-degree-of-freedom mathematical model for gimbaled 
proprotor (positive directions shown) . 
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WHIRL FLUTTER EQUATIONS OF MOTION FOR PROPROTORS WITH HINGED BLADES 


The linear equations of motion for the whirl flutter analysis of proprotors 
with hinged blades are given in reference 1 in the matrix form 


M(5> 

+ «[[b] + 

[D]]{4> . [i)2 [p] . 

M] * 

s 

II 



(B1) 

By arranging the four degrees of freedom indicated in 
order in the matrix the square matrices appearing 

appendix 5 of ref. 1) are given by 

figure B1 in sequential 
in equation (B1) (from 





Ip 

0 


I2 

0 




s 

II 


0 

l2 

I 

1 + *' ly 

0 

0 

I3 

l2 

0 


(B2) 



- 

0 



? 

0 

I3_ 





^ ^5 

HA3(ai 

- 32) 

A5 - eAij 

a 

iH(A3 

- eA2> 


[B] = K 

HA^Ca-i - a.2) 
A5 - eAi| 

a22H2A^ 

-a2H(A3 

+ A5 

- eA2) 

-a2H(A2 “ ^^2^ 
Ae 


A5 - 

eAij 

0 

(B3) 


a-fHCAg - 

eAg) 

A5 - 

eA^ 


0 



Ae 





-2I1 

0 

-212“ 







[D] = 

2I1 

CM 

cr 

0 

2I2 

0 





(B4) 



0 

-212 


-213 








2I2 

0 

2I3 
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[C] = K 


-aiH 3 Ai 


-h2ac 


-h 2(A3 - eA2> 


h2a^ 


-a 2 H 3 Ai 


a-|H(A3 - eA2) 
+ B3 tan 63 

(A5 - eA4) 

- a2HBi tan 63 


h 2(A3 - eA2) (B3 - eB2) tan 63 


-(A5 - eA4) 

+ aiHBi tan 63 

a2H(A3 - eA2) 

+ B3 tan 63 

-Ag 

(B3 - eB 2) tan 6^ 

(B 5)3 


[E] = 


0 

0 

0 

0 


Vt2(i, 


[S] = Q 2 


0 
0 
0 
0 

ai2R2M 

0 

0 

0 


0 

0 

eS 

0 


0 

0 

0 

eS 


Ip) 


V22(Ii + a22p2M 
0 
0 


ly) 


(B6) 


0 

V32I3 

Cr/J2 


0 

0 

-Cr/Q 

V32I3 


(B 7 ) 


The notation of reference 1 is maintained here for ease of comparison with 
the original work. However, several additional terms, not given in reference 1 , 
have been included here to represent the pylon Inertia (Ip and ly) and viscous- 
type structural damping for the pylon and rotor (Cq^, Cqg, and Cr). Also, to 


provide for a more realistic account of the lifting portion of the blade, the 
limits of integration in the aerodynamic integrals Ajj and Bjj were taken as 
TI1 and TI2 (as in the gimbaled proprotor equations) rather than e and 1 (as 
in ref. 1). 


3 lf flapping is defined about a line normal to the spanwise axis of the 
blade and lying in the plane of the hub as is customary, kinematic pitch-flap 
coupling appears as tan 63. Reference 1 obtained sin 63 because flapping 
was incorrectly defined there- as being about the skewed flapping hinge. 
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The viscous damping coefficients Cg, 
[s] are evaluated from the expressions; 


Cq^ , and Cqg appearing in and 


Cr = 21353^^1 + || 

Cq., = 2?pWp,o^P,o 
Cqg = 2 Cx<Dy,oIy,o 


(B8) 


where C3, Cp, and Cy damping ratios in flap, pitch, and yaw, respec- 
tively, wp o = 2TTfp o» ‘*>Y,o = 2TTfY o» ^P>o mass 

moments of inertia or the model in pitch and yaw with the blades locked to 
the forward flapping stops. 


Following reference 1 , introducing the nondimensional time T = Qt in 
equation (B1), and dividing by leads to 


w{i>" * [M * - rra * m 

where 



(B9) 


(BIO) 


By combining matrices in equation (B9), the equations of motion can be cast 
into the general form 


» Cc:{q}’ . M{q> = {0} 

where 

M = [A] 

[c] = [B] + [Dj i 

[K] = 05] ^ [E] +J_[S]| 


(BID 


(B12) 


By premultiplying equation (BID through by [k]”"* 

M -'M{q>" * M-' [C]{q}' - ra{q> = (o) (B13)" 


'^The equations of reference 1 employed here are cast into standard eigen- 
value form by first premultiplying the matrix equation by the inverse of [k] 
rather than by the inverse of [M] as in appendix A. The choice is arbitrary. 
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and introducing the four generalized velocities as auxiliary variables by means 
of the matrix identity 


{,}' - {q}’ = {0} 

equations (B13) and (B14) can be combined to yield 





[I] 

-W-i[c] 



(B14) 


(B15) 



By defining 

(x>= (B16) 

equation (B15) can be written in the compact form 

(x) = [a](x}’ (B17) 

where the definition of [a] follows from equation (B15). Substituting the 
assumed solution 


{x} = {x„)e>‘X 

into equation (B17) reduces equation (B17) to 
(Xo) = X[A]{Xo) 

or, by rearranging, to the familiar standard eigenvalue form 
[A]{Xo) □ X{Xo) 


where 


1/X 


(B18) 


(B19) 


(B20) 


(B21) 


Solution of equation (B20) leads to eight complex eigenvalues Xp and 
eigenvectors which occur in four complex conjugate pairs. The pth com- 


plex eigenvalue having a positive imaginary part has the form 

Xp = Up + iVp 

and can be interpreted as 


_ _ t*JpCp ^ ^ ‘^pX/'* ~ ^p^ 


^2 


a 


(B22) 


(B23) 


where cop is the undamped coupled frequency and Cp is the damping as a frac- 
tion of critical damping. 
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The complex vector associated with the eigenvalue Xp has the form 





{B2H) 


where the lower four elements in each modal vector of eight elements 

define the mode shape. As before, the relative amplitude and phasing between 
the pylon pitching and yawing motions in a given mode can be ascertained by 
converting the complex elements of the mode to polar form and then normalizing 
on one of them. 
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Degrees of freedom 

Pylon pitch rotation 
Q2 Pylon yaw rotation 

qg Longitudinal flapping of proprotor disc 
q^ Lateral flapping of proprotor disc 



Figure B1 Four— degree—of— freedom mathematical model for proprotor with 
offset flapping hinges (positive directions shown) . 
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APPENDIX C 


DETERMINATION OF EQUIVALENT HUB SPRING 
Blade Flapping Natural Frequency 


Gimbaled proprotor .- An expression for the blade flapping natural frequency 
can be derived from the equations of motion developed in appendix A by consid- 
ering the degenerate case obtained by retaining only the two cyclic flapping 
degrees of freedom of the proprotor disc. If precone and structural damping 
are neglected, the resulting equations for the case of a symmetric hub restraint 
reduce to 


Ir® 1 2Ipftbi + K{ja-| 


= I - 


33a 1 


IpBl - 2Ipflai + K{]b-| 



tan 63 


tan 63 



(Cla) 


(Clb) 


Now, from reference 15, the flapping angle B of the reference blade is related 
to the cyclic flapping degrees of freedom of the proprotor disc by 


3 = a-| sin \1> - b-| cos \}i 


(C2) 


where \j> = Jit is the azimuth angle of the reference blade. Multiplying equa- 
tion (Cla) by sin and equation (Clb) by cos subtracting the second equa- 
tion from the first, and making use of equation (C2) and the additional relations 

3 = a-) sin - b-| cos + a-|Si cos + biJ2 sin \|» 

\ (C3) 

3 = a-| sin - b-| cos \|> + 2aiJi cos + 26-|J2 sin ? 

- aiJ)2 sin + biJ)2 cos \|) ^ 


which follow from equation (C2), the resulting equation can be put into the 
form 


3+1 Y^A53 + 



+ ^^ + I YJi^B3 tan 63)3 = 

Ir 2 


0 


(C4) 


From equation (C4) the undamped flapping natural frequency 0)g, nondimension- 
alized by the rotor rotational speed, is given by 


(03 



+ 3 l 

lRft2 


^ YB3 tan 63^ 


1/2 


(C5) 
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and the aerodynamic flap damping by 


= 


YA5 

4Q3 


(C6) 


In the absence of hub restraint and pitch-flap coupling, the blade flapping 
natural frequency in the rotating system is equal to the rotational speed 
that is, one per revolution. The use of hub restraint in combination with posi- 
tive or negative 63 raises or lowers the natural frequency from one per revo- 
lution. The effect of 63 on the blade flapping frequency is seen to be that of 
an aerodynamic spring which increases or decreases the frequency depending on 
whether 63 is positive or negative. 


Proprotor with offset flapping hinges .- By using the equations of motion 
given in appendix B for a proprotor having hinged blades , the expressions anal- 
ogous to equations (C 5 ) and (C6) can be shown to be 


(03 = 


/l + + 

I ^3 


V 3 ^ 


^(B 3 - 


eBg) tan 63 


1/2 


(C7) 


C3 = ^(A5 - 2 eA 4 + (C8) 

where 

Ye = PacR^ (C9) 

2I3/N 


is the Lock number defined by using the mass moment of inertia of the blade 
about the offset flapping hinge rather than about the center line of the shaft. 
By comparing equations (C 5 ) and (C6) with equations (C 7 ) and (C8), it is seen 
that several additional terms appear for the proprotor having noncentrally 
hinged blades. In the absence of flapping restraint, hinge offset, and pitch- 
flap coupling, the flapping natural frequency is one per revolution. The use 
of flapping restraint or hinge offset in combination with positive or negative 
63 raises or lowers the natural frequency from one per revolution. 


Basis of Establishing Equivalence 

A gimbaled proprotor is essentially a proprotor having centrally hinged 
blades. On the supposition that an "equivalent” hub spring can be used to rep- 
resent the inertial gind aerodynamic restoring moment effects associated with an 
offset flapping hinge, analyses developed for the gimbaled proprotor would then 
also be applicable to proprotors having blades with offset flapping hinges, at 
least for preliminary design calculations. An equivalent hub spring can be 
established by requiring that the flapping natural frequency be preserved, that 
is, that 

^3! gimbaled = ^Blhinged (CIO) 
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This requirement implies that 


% 

lRft2 


+ J. Y ^3 tan 63 = + .^£(63 - € 82 ) tan 63 


By approximating Ye ty Y> equation (C11) leads to 

= ®§ + Vo 2 _ 1 Y 6 B 2 tan So 
lRft2 I 3 3 2 ' 2 3 


(C11) 


(C12) 


so that the equivalent hub spring is given by 

Kfl.eq = + '’3^ - ^ Y£B 2 tan (C13) 

For a given flapping hinge proprotor configuration, and a selected rotational 
speed and airspeed, the equivalent hub restraint is easily determined from equa- 
tion (CI 3 ). It should be noted that a new equivalent hub restraint must be 
determined for each ft,V combination. 


The equivalent hub restraint given by equation (Cl 3) was established by 
requiring that the flapping natural frequency be preserved while retaining both 
the dynamic and aerodynamic terms in the expressions for (5g. An equivalent 
hub spring established on the basis of preserving the in-vacuum flapping natural 
frequency by discarding the aerodynamic terms in equation (C11) would not be 
satisfactory. 
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TABLE I.- PHYSICAL PROPERTIES OF MODEL 



SI Units 


U.S. Customary Units 


e/R = 0.05 

e/R □ 0.13 

e/R = 0.05 

e/R = 0.13 

Rotor: a 

Blade radius 

0.744 m 

0.805 m 

2.44 ft 

2.64 ft 

Blade chord 

0.0902 m 

0.0902 m 

0.296 ft 

0.296 ft 

Blade mass 

0.533 kg 

0.508 kg 

0.0365 slug 

0.0348 slug 

Blade static moment about 
flapping hinge 

0.111 kg-m 

0.106 kg-m 

0.025 slug-ft 

0.0238 slug-ft 

Blade moment of inertia 
about flapping hinge . . . 

0.0439 kg-m2 

0.0436 kg-m^ 

0.0324 slug-ft2 

0.0322 slug-ft2 

Distance from pylon pitch 
axis to rotor hub .... 

0.320 m 

0.320 m 

1.05 ft 

1.05 ft 

Distance from pylon yaw 
axis to rotor hub .... 

0.277 m 

0.277 m 

0.91 ft 

0.91 ft 

Pylon: 

Mass effective in 

pitch 

3.37 (3.46) kg 

3.87 kg 

0.231 (0.237) slug 

0.265 slug 

Mass effective in yaw . . . ■ 

3.01 kg 

3.08 kg 

0.206 slug 

0.211 slug 

Pitch inertia about center 
of gravity 

0.0496 (0.0526) kg-rm2 

0.0565 kg-m^ 

0.0366 (0.0388) slug-ft2 

0.0417 slug-ft2 

Yaw inertia about center of 
gravity 

0.0343 kg-m2 

0.0502 kg-m^ 

0.0253 slug-ft2 

0.0370 slug-ft2 

Distance from pylon pitch 
axis to pylon center of 
gravity 

0.212 (0.207) m 

0.187 m 

0.694 (0.679) ft 

0.614 ft 

Distance from pylon yaw 
axis to pylon center of 
gravity 

0.193 m 

0.168 m 

0.632 ft 

0.552 ft 


®Blade inertial properties varied slightly from blade to blade. Values shovm are average for the three 
blades . 

^Pylon inertial properties include rotor hub but exclude blades. 

^Values in parentheses reflect changes to pitch spring mechanism for runs 62 to 68. 



TABLE II.- EXPERIMENTAL WHIRL FLUTTER RESULTS 


Run 

Point 

e/R 

— 

«3, 

deg 

Hz 

fY,o» 

Hz 

2Cp 

2 Cy 

n, 

Hz 

o)p/Q 


Vp/«R 

top /ft 

Whirl 

+ or - 

Ap 

4>p, 

deg 

40 

26 

0.05 

20 

5.60 

5.80 

0.024 

0.051 

i 13.3 

0.444 

0.462 

0.78 

0.42 

+ 

1.66 

72 


27 

.05 

20 

5.60 

i 5.80 

.024 

.051 

14.0 

.422 

.439 

.75 

.39 

+ 

1.47 

82 

41 

i 6 

.05 

20 

5.40 

5.50 

.024 

.051 

8.0 

.712 

.729 

1.22 

.66 


1.48 

61 


9 

.05 

20 

5.40 

5.50 

.024 

.051 

11.1 

.513 

.525 

.88 

.50 

+ 

1.64 

65 


13 

.05 

20 

5.40 

5.50 

.024 

.051 

17.0 

.335 

.343 

.60 

.34 

+ 

1.54 

72 

42 

6 

.05 

20 

, 3.76 

3.82 

1 .012 

! 

.048 

13.9 

.285 

.292 

.56 

.29 

+ 

2.11 

73 


8 

.05 

20 

3.76 

3.82 

.012 

.048 

8.0 

.496 

.507 

.76 

.50 


1.24 

90 


12 

.05 

20 

3.76 

3.82 

.012 

.048 

: 11.0 

.361 

.368 

.60 

.36 

+ 

1.27 

76 


15 

.05 

20 

3.76 

3.82 

.012 

.048 

17.0 

.233 

, .238 

.44 

.24 


1.78 

65 

43 

5 

.05 

20 

3.74 

5.70 

.012 

.046 

8.0 

.496 

.760 

1.12 

.24 


.58 

-54 


7 

.05 

20 

3.74 

5.70 

.012 

.046 

11.0 

.360 

; .551 

.81 

.24 

- 

.35 

-99 


8 

.05 

20 

3.74 

5.70 

: .012 

.046 

14.0 

, .281 

: .430 

1 .56 

.23 

- 

.19 

-72 


12 

.05 

20 

3.74 

5.70 

i 

.012 

.046 

16.8 

.235 

.359 

.47 

.20 

- 

.24 

-60 

44 

3 

.05 

30 

3.74 

I 

1 3.90 

1 .012 

.048 

8.0 

.490 

.513 

.56 

.56 

+ 

1.36 

65 


4 

.05 

30 

3.74 

1 3.90 

j .012 

.048 

J 

10.8 

.365 

' .383 

.40 

.41 

+ 

1.76 

61 

45 

4 

.05 

30 

5.64 

5.54 

’ .019 

1 

.041 

8.0 

. .748 

.739 

.78 

' .77 

+ 

1.36 

87 


8 

.05 

30 

5.64 

5.54 

.019 

.041 

, 10.9 

.545 

. .538 

.55 

.56 

•f 

1.34 

86 


10 

.05 

30 

5.64 

5.54 

' .019 

.041 

13.8 

.431 

.425 

.41 

.46 

+ 

1.43 

81 


11 

.05 

30 

5.64 

5.54 

.019 

.041 

15.8 

1 

.377 

.372 

.38 

.40 

+ 

1.32 

94 

46 

7 

. .05 

30 

6.76 

6.78 

.018 

.021 

1 

' 10.9 

.654 

.659 

.64 

.70 


2.05 

68 


14 

.05 

1 30 

6.76 

6.78 

.018 

.021 

14.0 

.511 

' .515 

.50 

.54 

® + 

1.96 

68 


15 

, -05 

i 30 

6.76 

6.78 

, .018 

.021 

17.0 

.421 

.424 

.40 

.46 

+ 

1.97 

97 


18 

: .05 

30 

6.76 

6,78 

.018 

.021 

10.8 

.660 

.666 

1 .66 

.70 

+ 

1.95 

75 

' 48 

9 

.05 

' 6.75 

6.76 

6.78 

.013 

.037 

10.9 

.363 

.361 

, .82 

.28 

+ 

1.02 

83 


11 

: .05 

i 6.75 

6.76 

6.78 

.013 

.037 

13.8 

.286 

.286 

1 -59 

.23 

- 

1.24 

-61 

' 

I 

.05 

6.75 

6.76 

j 

6.78 

-1 

.013 

.037 

1 

16.9 

J 

.235 

1 

.234 

i .51 

J 

.18 

-1 

-i 

1 2.06 

-69 



TABLE II.- Continued 


UJ 


Run 

Point 

e/R 

« 3 , 

deg 

fp,o> 

Hz 

fY,o» 

Hz 

1 

2 ^p 

2 Cy 

n, 

Hz 

(A)p/fi 


Vp/JiR 

a)p/il 

Whirl 

+ or - 

Ap 

<t>F> 

deg 

50 

11 

0.05 

6.75 

5.54 

5.34 

0.005 

0.049 

16.9 

0.347 

0.336 

0.74 

0.28 

+ 

1 .07 

80 

51 

14 

.13 

20 

5.30 

4.65 

.011 

.055 

10.8 

.518 

.456 

.92 

.21 

_ 

NA 

NA 


22 

.13 

20 

5.30 

4.65 

.011 

.055 

14.1 

.397 

.350 

.72 

.16 

- 

NA 

NA . 

52 

9 

.13 

20 

5.30 

4.65 

.027 

.088 

16.6 

.335 

.296 

.65 

.13 

- 

2.01 

-76 

53 

9 

.13 

20 

3.72 

3.56 

.009 

.034 

10.9 

.360 

.347 

.73 

.16 


NA 

NA 


12 

.13 

20 

3.72 

3.56 

.009 

.034 

13.9 

.283 

.271 

.60 

.12 


NA 

NA 


17 

.13 

20 

3.72 

3.56 

.009 

.034 

16.8 

.234 

.225 

.50 

.10 

- 

NA 

NA 


19 

.13 

20 

3.72 

3.56 

.009 

.034 

7.8 

.506 

.486 

.92 

.22 

- 

NA 

NA 

54 

7 

.13 

10.5 

3.72 

3.52 

.013 

.040 

8.0 

.410 

.466 

1.01 

.19 

- 

NA 

NA 

55 

13 

.05 

20 

3.70 

3.58 

.013 

.040 

11.1 

.351 

.341 

.42 

.38 

+ 

.68 

98 


15 

.05 

20 

3.70 

3.58 

.013 

.040 

13.8 

.283 

.274 

.34 

.30 

+ 

.50 

81 

56 

2 

.05 

20 

3.75 

3.56 

.013 

.037 

10.9 

.363 

.347 

.38 

.38 

+ 

.66 

87 


3 

.05 

20 

3.75 

3.56 

.013 

.037 

13.2 

.300 

.286 

.34 

.31 

+ 

.96 

75 


7 

.05 

20 

3.75 

3.56 

.013 

.037 

8.0 

.495 

.472 

.53 

.51 

+ 

.65 

101 

57 

14 

.05 

20 

3.72 

3.82 

.008 

.031 

8.0 

.490 

.507 

.68 

.50 

+ 

.75 

64 


15 

.05 

20 

3.72 

3.82 

.008 

.031 

10.1 

.353 

.365 

.44 

.37 

+ 

.82 

63 


17 

.05 

20 

3.72 

3.82 

.008 

.031 

14.2 

.276 

.285 

. 34 + 

.30 

+ 

.73 

63 


18 

.05 

20 

3.72 

3.82 

.008 

.031 

15.2 

.257 

.265 

. 34 - 

.27 

+ 

.67 

95 

58 

8 

.05 

20 

3.72 

5.60 

.008 

.061 

8.0 

.487 

.727 

1.08 

.23 

_ 

.28 

-87 


12 

.05 

20 

3.72 

5.60 

.008 

.061 

11.1 

.353 

.527 

.79 

.23 

- 

.19 

-88 


15 

.05 

20 

3.72 

5.60 

.008 

.061 

14.2 

.276 

.412 

.55 

.22 

- 

.15 

-60 

60 

4 

.05 

20 

5.48 

5.44 

.010 

.049 

11.0 

.525 

.524 

.61 

.54 

+ 

.46 

86 


6 

.05 

20 

5.48 

5.44 

.010 

.049 

14.0 

.412 

.411 

.46 

.42 

+ 

.45 

97 


8 

.05 

20 

5.48 

5.44 

.010 

.049 

16.9 

.342 

.341 

.36 

.37 

+ 

.53 

75 


TABLE II.- Concluded 


U) 

00 


Run 

Point 

e/R 

63, 

deg 

Hz 

^■7,0* 

Hz 

2Cp 

2Cy 

n, 

Hz 

(A)p/{2 

a>y/fl 

Vp/QR 

(Op/^2 

Whirl 

+ or - 

Ap 

♦p, 

deg 

62 

5 

0.05 

20 

5.40 

5.68 

0.017 

0.034 

8.1 

0.704 

0.743 

1.13 

0.68 

+ 

NA 

NA 


9 

.05 ! 

20 

5.40 

5.68 

.017 

.034 

11.2 

.508 

.537 ' 

.76 ! 

.51 

+ 

NA 

NA 


10 

.05 

20 

5.40 

5.68 

.017 

.034 

14.1 

.404 

.427 

.56 

.42 

+ 

NA 

NA 


11 

.05 

20 

5.40 

5.68 

.017 

.034 

17.3 

.329 

.348 

.46 

.35 

+ 

NA J 

NA 

i 63 

12 

.05 

20 

7.06 

7.40 

.022 

.034 

11.2 

.664 

.700 

1 1.07 

.64 

+ 

.52 

83 

1 

13 

.05 

20 

7.06 

7.40 

,022 

.034 

14.2 

.526 

.554 

.81 

.50 


.67 

65 


20 

.05 

20 

7.06 

7.40 

.022 

.034 

17.4 

.429 

.452 

.66 

.44 

+ 

.63 

90 


27 

.05 

20 

7.06 

7.40 

,022 

.034 

9.0 

.827 

.871 

1.36 

.78 

I + 

.58 ! 65 

64 

6 

.05 

20 

4.00 

3.96 

.024 

.044 

: 8.2 

.517 

.514 

.67 

.53 

+ 

.61 

104 


8 

.05 

20 

I 4.00 

3.96 

.024 

.044 

11.2 

.377 

.375 

.46 

.40 

+ 

.62 

65 

1 

9 

.05 

20 

4.00 

3.96 

.024 

.044 

14.2 

.298 

.297 

j .36 

.32 

+ 

.60 

72 


11 

.05 

20 

4.00 

3.96 

.024 

.044 

15.2 

.278 

.277 

.35 

.31 

.+ 

.60 

79 

' 65 

i 

1 15 

.05 

20 

i 4.60 

3.96 

.024 

1 

i .044 

i 8.2 

.592 

.512 

\ 

1.06 

j .50 

+ 

NA 

NA 

' 

i 22 

' .05 

1 20 

4.60 

3.96 

, .024 

1 .044 

' 11.1 

.437 

.378 

' .75 

! .36 

+ 

NA 

NA 


31 

.05 

i 20 

4.60 

3.96 

.024 

i .044 

i 14.3 

i .341 

.294 

.61 

! .29 

+ 

NA 

NA 


' 36 

.05 

i 20 

, 4.60 

3.96 

1 .024 

; .044 

1 17.2 

.282 

.244 

, .50 
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Figure 4 .- Hub and control system geometry. Figure is not to scale 
and pitch horn is in plane of hub. 63 □ tan~^ (d-|/d2). 
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Figure 8.- Summary plot of all experimentally identified whirl flutter points. 


50 


Flutter inflow ratio, Vjj/flR Flutter inflow ratio, 


I 


Theory E3q)eriment 

^ Reference 15 O Forward whirl 

~ — — Reference 1 # Backward whirl 



0 .2 .4 .6 .8 0 .2 .4 .6 .8 

Pylon frequency ratio, Pylon frequency ratio, 



0 .2 .4 .6 .8 0 .2 .4 .6 .8 

Pylon frequency ratio, Pylon frequency ratio, o}p/Sl 


Figure 9.- Run-by-run correlation of experimental and theoretical 

flutter inflow ratios. 
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Figure 9.- Continued. 
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Figure 9.- Continued. 
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Continued. 
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Figure 9.- Continued. 




Flutter inflow ratio, Flutter inflow ratio, 


Theory 


Reference 15 
Reference 1 



Experiment 

O Forward whirl 
# Backward whirl 





Figure 9.- Concluded. 
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Figure 10.- Run-by-run correlation of experimental and 
theoretical flutter frequencies. 
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■Figure 10.- Continued. 


58 



Flutter frequency ratio, Flutter frequency ratio, 


Theory 


Experiment 

Reference 15 O Forward whirl 

Reference 1 ® Backward whirl 



0 .2 .4 .6 .8 

Pylon frequency ratio, 



Pylon frequency ratio, o^u/ 



0 .2 .4 .6 .8 


Pylon frequency ratio, 



0 .2 .4 .6 .8 

Pylon frequency ratio, 


Figure 10.- Continued. 
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Figure 10. 
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Figure 10.- Concluded. 
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Figure 11.- Comparison of measured pylon yaw-to-pitch amplitude ratio 
at flutter with theory of reference 1. 
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Figure 12.- Comparison of measured pylon yaw-to-pitch phase angle 
at flutter with theory of reference 1 . 
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Figure 14.- Effect of hinge offset on whirl flutter. 
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Figure 15.- Effect of pylon support stiffness on whirl flutter. For toy < t»)p, 
flutter inflow ratio is plotted against 
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Figure 16.- Effect of flap hinge damping on forward whirl flutter 
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Figure 17.- Effect of pylon damping on forward whirl flutter. 
Symmetric pylon frequency . 
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Figure 19.- Effect of pylon damping on backward whirl flutter. 
Nonsymmetric pylon frequency. 
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Figure 20.- Effect of flap hinge damping on backward whirl flutter 

Symmetric pylon frequency. 
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